Expansion of a Fermi Cloud in the BCS-BEC Crossover 
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We study the free expansion of a dilute two-component Fermi gas with attractive interspecies 
interaction in the BCS-BEC crossover. We apply a time-dependent parameter-free density-functional 
theory by using two choices of the equation of state: an analytic formula based on Monte Carlo data 
and the mean-field equation of state resulting from the extended BCS equations. The calculated 
axial and transverse radii and the aspect ratio of the expanding cloud are compared to experimental 
data on vapors of 6 Li atoms. Remarkably, the mean-field theory shows a better agreement with the 
experiments than the theory based on the Monte Carlo equation of state. Both theories predict a 
measurable dependence of the aspect ratio on expansion time and on scattering length. 



PACS numbers: PACS Numbers: 03.75.Kk 

I. INTRODUCTION 

Current experiments with cold vapors of 6 Li and 40 K 
atoms can operate in the regime of deep Fermi degener- 
acy. The available experimental data on two-hyperfine- 
component Fermi gases are concentrated across a Fesh- 
bach resonance, where the s-wave scattering length ojr 
of the interatomic potential varies from large negative to 
large positive values Q, 0, EJ and where a crossover 
from a Bardeen-Cooper-Schrieffer (BCS) superfluid to a 
Bose-Einstein condensate (BEC) of molecular pairs has 
been predicted |a,|(|l3- In these experiments, the Fermi 
cloud is dilute because the effective range Rq of the inter- 
action is much smaller than the mean interparticle dis- 
tance, i.e. kpRo *C 1 where kp = (Stt 2 ?!) 1 / 3 is the Fermi 
wave vector and n is the gas number density. Even in 
this dilute regime the s-wavc scattering length ap can 
be made very large: the interaction parameter kpap di- 
verges and changes sign at a Feshbach resonance, despite 
kpR remaining small 0,0 IS II]- 

Recent experimental and theoretical investigations 
studied the densit y p rofil es 1^, IToj. collective exci- 
tations [IJllElM Q lial3nm_condensate frac- 
tion pLlral2^n2lfand vortices |22Tl23| of the fermion 
cloud through the BCS-BEC crossover. In this letter 
we analyze the free expansion of the Fermi gas through 
this crossover by using a parameter- free time-dependent 
density- functional theory jlalr7| | based on the bulk equa- 
tion of state of the superfluid, and including a quantum- 
pressure term. We adopt two possible equations of state: 
a reliable analytical interpolating formula based on bulk 
Monte Carlo results an d the mean-field equation of 
state based on extendend BCS equations [EEISE!- Ex- 
perimentally, the free expansion of superfluid 6 Li clouds 
was observed by O'Hara et al. 0] and by Bourdel et al. 
4] . The comparison of our theory with these experimen- 
tal data shows that the effects of interaction could be 
detected during the expansion if the thermal component 
was negligible. In addition, by using local scaling equa- 
tions, we investigate the long-time dynamics of the Fermi 
gas predicting novel and measurable effects of interaction 
on the time evolution of the expansion process. 



II. THEORY 

To describe the dynamics of a zero-temperature Fermi 
cloud in the external potential U(r,t) we use a hydro- 
dynamic model with a von Weizsacker quantum-pressure 
term. This theoretical approach is expected to be reli- 
able for studying the collective dynamics of the Fermi 
gas 0,^3- The action functional A[ijj} of the theory de- 
pends on the superfluid order parameter ip(r, t) as follows 

A = J dt d 3 r £(ip,dtip,Vijj) , (1) 
where the Lagrangian density reads 

c = in rdt^P + |-V*vV - u\ip\ 2 - sm 2 M 2 - (2) 

2m 

£ represents the bulk energy per particle of the system, 
which is conveniently expressed as a function of the num- 
ber density n = |^/)| 2 by the following equation: 

E(n) = \ *J± f(y) , (3) 
5 2m 

where f(y) is a function of the inverse interaction pa- 
rameter y = (kpap) -1 ■ In the weakly attractive regime 
(y <C —I) one expects a BCS Fermi gas of weakly bound 
Cooper pairs where the superfluid gap energy A is expo- 
nentially small. In the so-called unitarity limit (y = 0) 
one expects that the energy per particle is proportional to 
that of a non-interacting Fermi gas with a n-independent 
coefficient /(0) = 0.42 24]. In the weak-coupling BEC 
regime (y I), a weakly repulsive Bose gas of dimcrs 
of mass rnj = 2m and density rig = n/2 is expected. 
Such Bose-condensed molecules interact with a positive 
scattering length as = 0.6a_p [2^, [2^] . The function f(y) 
is modelled by the analytical formula 

f(y) = ai-a 2 arctan ( a 3 y ^ + pr ) (4) 

V P2 + \y\J 

recently derived from Monte Carlo (MC) simulations 
[2^,|23 and asymptotic expressions. Table I of Ref. 
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reports the values of the interpolating a\, 0,2, 123, and 
/3 2 . 

In the present investigation we take an axially symm- 
metric harmonic potential as confining trap 

U(r, t) = ^ [w„(t) V + V 2 ) + ^{tfz 2 } , (5) 

where tOj(t) — ujjQ{—t), with j = 1,2,3 = p,p,z and 
O(i) is the step function, so that, after the external trap 
is switched off at t > 0, the Fermi cloud performs a free 
expansion. The Euler-Lagrange equation for the field 
"0(r, t) is obtained by minimizing the action functional 
of Eqs. (1,2). This leads to a time-dependent nonlinear 
Schrodinger equation (TDNLSE): 



2m 



V 2 + [/ + M (|V| 2 ) 



(6) 



The nonlinear term fi is the bulk chemical potential of 
the system. Like the energy £ of Eq. (3), also the bulk 
chemical potential /1 is a function of the number den- 
sity n. The MC chemical potential is related to the MC 
energy by the thermodynamical formula 
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Instead of the MC equation of state (7) based on (4), 
in the TDNLSE (6) one can plug the mean-field equa- 
tion of state, obtained from the extendend BCS (EBCS) 
equations I n this scheme, the chemical poten- 

tial \x and the gap energy A of the uniform Fermi gas are 
found by solving the following EBCS equations 
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where I\ (x) and I2 (x) are two monotonic functions which 
can be expressed in terms of elliptic integrals [2(], . By 
solving these two EBCS equations we obtain the chemi- 
cal potential \i as a function of n and a^, which can be 
inserted into the TDNLSE (6). 



III. EXPANSION OF FERMI GAS AND 
SCALING EQUATIONS 

The free expansion of a droplet of 1.5- 10 5 6 Li atoms in 
the unitarity limit (y ~ 0) was investigated experimen- 
tally in Ref. pj. The harmonic potential is anisotropic 
with A = tjJ z jtjJp = 0.035. The scattering length for the 
applied magnetic field B = 910 G is op = —0.38 pm 
= —0.14 a z [33, which correspond to y = —0.16 at 
the droplet center. Here a z — [/i/fmuj] 1 ' 2 . Fig- 
ure 1(a) compares the observed full width half maximum 
(FWHM) of the transverse and axial size of the expand- 
ing cloud as a function of time with the ones obtained 
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FIG. 1: (Color online). Expansion of a cloud of 1.5 • 10 5 
6 Li atoms released from a trap as realized in Ref. Q, with 
anisotropy A = lo z /lo p = 0.035. (a) Transverse and axial 
radii of the 6 Li atomic cloud close to the unitarity limit: 
hf = — 0.14a z which corresponds to y = —0.16. (b) As- 
pect ratio of the cloud as function of the time t. Circles 
and squares: experimental data of Ref. [lj; dot-dashed lines: 
TDNLSE with the MC equation of state; solid lines: LSE 
with the MC equation of state; dashed lines: LSE with the 
EBCS equation of state. 



by numerical integration of the TDNLSE, based on the 
MC and EBCS equation of state. The TDNLSE is solved 
numerically by using a finite-difference algorithm §3l| on 
a real-space grid. 

Figure 1(a) shows that the expanding gas accelerates 
more strongly in the radial direction, where the confine- 
ment is tighter, than axially. Accordingly, the cloud un- 
dergoes a shape transition: from a cigar to disk. This 
is a consequence of superfluidity and interaction: a non- 
interacting or a normal Fermi gas would undergo a ballis- 
tic expansion, leading eventually to a spherical shape . 
In Fig. 1 the TDNLSE results are plotted as dot-dashed 
lines and show a fair agreement with the experimental 
data. It is important to observe that the present theory 
does not rely on any fitting parameters, while the model 
curves shown in Ref. [lj critically depend on the choice 
of the initial widths. The initial profile for the TDNLSE 
is obtained by running the code integrating Eq. (6) in 
imaginary time until the confined ground-state is filtered 
out. Figure 1(b) plots the droplet aspect ratio showing 
that the theory overestimates the experimental data. 

From the TDNLSE one can deduce Landau's hy- 
drodynamic equations of superfluids at zero tempera- 
ture, by setting ip{r,t) = \/n(r, i)e lS(r '*) and v(r, t) = 
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(H/m)'VS(r, t), and neglecting the quantum-pressure 
term (— S 2 V 2 v / n)/(2m v / n), that is expected to be com- 
parably small for a large number N of particles 0, 0, 
|29j. These hydrodynamic equations are 

d t n + V • (nv) = , (10) 



m <9 t v + V ( /i(n) + U (r, <) + —mv 2 



0. (11) 



In this approximation, the stationary state in the 
trap is given by the Thomas- Fermi profile ?io( r ) = 
(p, — U(r, 0)). Here p,, the chemical potential of 
the inhomogeneous system, is fixed by the normaliza- 
tion condition N = J d 3 r nn(r). We impose that 
the hydrodynamic equations satisfy the scaling solutions 
n(r,f) = n (x/b 1 (t),y/b 2 (t),z/b 3 (t))/b{t) and v(r,t) = 

(xb 1 (t)/bi(t),yb 2 (t)/b 2 (t),zb 3 (j:)/b3(t)y where 6(t) = 

rifc=i &fcC0- We obtain three differential equations for the 
scaling variables bj(t), with j = 1,2, 3 = p,p,z. These 
scaling differential equations depend also on the space 
vector r. Only if the chemical potential satisfies a poly- 
tropic power law u(n) = Cn 1 then the space dependence 
drops out |14i Il6|. In our problem /x(n) is not a power 
law but we expect that the dynamics can be well approx- 
imated by evaluating the scaling differential equations at 
the center (r = 0) of the cloud j^]. In this case the 
variables bj (t) satisfy the local scaling equations (LSE) 
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The coupled ordinary differential equations (12) are in- 
tegrated accurately and efficiently to arbitrary time by 
standard algorithms. We check the reliability of the LSE 
approach by comparing their numerical solutions to the 
expansion obtained by using the full TDNLSE (6), both 
within the MC equation of state (7). Figure 1 reports the 
LSE results as solid lines, clearly showing that the LSE 
are extremely accurate: solid lines are practically super- 
imposed to dot-dashed lines (relative difference -C 1%, 
see inset of Fig. 1(b)). Figure 1(a) also reports the trans- 
verse and axial radii obtained by solving the LSE with 
the chemical potential /i(n) given by the EBCS equations 
(8-9). Remarkably the mean- field EBCS results are closer 
to the experimental data than the MC results. The two 
theories essentially coincide for the aspect ratio. 

In Ref. Q the free expansion of 7 • 10 4 cold 6 Li atoms 
was studied for different values of y = (kpap)^ 1 around 
the Feshbach resonance (y = 0). Unfortunately, in this 
experiment the thermal component is not negligible and 
thus the comparison with the present T = theory is not 
fully satisfactory. Figure 2 compares the experimental 
data of Ref. [j] with the LSE based on both MC and 
EBCS equation of state. Figure 2 shows that the aspect 
ratio predicted by the two T = theories exceeds the 
finite-temperature experimental results. This is not sur- 
prising because the thermal component tends to suppress 
the hydrodynamic expansion of the superflud. On the 



1.8 
1.6 
1.4 
1.2 



T 



t= 1.4 ms 



(a) 



/ / □ 

MC eqn of state rfKti D n 




^ [] T. Bourdel et al. 
b 




EBCS eqn of state 



(b) 



FIG. 2: (Color online). Properties of a 6 Li cloud after 1.4 ms 
expansion from the trap realized in Ref. Q, of anisotrpy A = 
lu z /lu p — 0.34. (a) Aspect ratio of the 6 Li atomic cloud as a 
function of the inverse interaction parameter y — (/s^of) -1 ■ 
(b) Released energy of the same cloud defined as in Ref. Q 
based on the rms widths of the cloud, (c) Actual released 
energy of the atomic cloud. Squares report the experimental 
data of Ref. 3; solid lines: LSE with the MC equation of 
state; dashed lines: LSE with the EBCS equation of state. 



other hand, the released energy of the atomic gas is well 
described by the two T = theories, and again the mean- 
field theory seems more accurate, also probably due to 
the thermal component. For completeness, Fig. 2(c) re- 
ports the actual released energy j_d 3 r no(r) £ (no(r)). 

In the two experiments of Ref. [jj and Ref. [j] the time 
evolution is sufficiently short for a full TDNLSE simu- 
lation. It would be computationally impractical to in- 
tegrate the TDNLSE for times much longer than uj^ 1 , 
where u>h = (wpW 2 )'' 3 - As the LSE are very reliable at 
small and intermediate times, we use them to investigate 
the time evolution of the Fermi cloud for longer times. 
Figure 3 shows the aspect ratio of the expanding cloud as 
a function of the inverse interaction parameter y at sub- 
sequent time intervals. At t = the aspect ratio equals 
the trap anisotropy A = 0.34. According to these calcula- 
tions, during the cloud expansion the aspect ratio in the 
BCS regime (y <C —1) is measurably different from the 
one of the BEC regime (y ^> 1). Thus the free expansion 
enables one to recognize the regime involved. Figure 3 
predicts a novel interesting effect: initially (tuiu ^5 3) 
the cloud aspect ratio evolves faster in the BCS region, 
but then at some intermediate time (here tuin — 4) the 
BEC side reaches and eventually overtakes the BCS side 
at larger times (here twjj > 5). Of course, the detailed 
sequence of deformations depends on the experimental 
conditions and in particular on the initial anisotropy, but 
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FIG. 3: (Color online). Successive frames of the aspect ra- 
tio of the Fermi gas as a function of the inverse interaction 
parameter y = (^fcsf) -1 in the experimental conditions of 
Ref. Q- At t = the Fermi cloud is cigar-shaped with 
a constant aspect ratio equal to the initial trap anisotropy 
A = ui z /ujp — 0.34. Solid lines: LSE with the MC equation of 
state; dashed lines: LSE with the EBCS equation of state. 



for the expansion of any initially cigar-shaped interacting 
fermionic cloud. Similarly, starting from a disk-shaped 
cloud (A > 1), the aspect ratio reduces more quickly ini- 
tially on the BCS side, and later on the BEC side. 



IV. DISCUSSION 

Comparison of the EBCS (dashed lines) and MC (solid 
lines) data shows that beyond mean-field effects do not 
alter qualitatively the general trend, but they affect the 
aspect ratio quantitatively, to an extent which could be 
appreciated by very accurate experiments carried out at 
extremely low temperature. In particular, the mean- field 
curves flatten to the asymptotic values (for \y\ ^> 1) 
closer to the unitary limit than the MC ones. Present-day 
experimental data, including measurements of collective 
oscillation frequencies HEflHElElEEIll, are 
equally well compatible with the EBCS mean field and 
the MC-based analysis accounting for beyond mean-field 
effects. New experiments could shed light on these cor- 
relation effects and verify the predictions of the present 
calculations. 
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